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It has recently been suggested that an exactly solvable problem characterized by a new quantum 
number may underlie the electronic shell structure observed in the mass spectra ol medium-sized 
sodium clusters. We investigate whether the conjectured quantum number 3n + 1 bears a similarity 
to the quantum numbers n + l and 2n + l, which characterize the hydrogen problem and the isotropic 
harmonic oscillator in three dimensions. 
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I. INTRODUCTION 



The electronic shell structure observed in mass spec- 
tra of metal clusters can be rather well described by an 
independent electron model. Treating the valence elec- 
trons of the cluster as an ideal Fermi gas moving in an 
external potential V(r), a shell structure arises from the 
discreteness of energy levels for particles confined to a fi- 
nite volume. This approach forms the basis for the semi- 
classical description of the electronic shell and supershell 
structureErul Suitable model potentials can be obtained 
from self-consistent density functional calculationsJHI It 
can be verified that the electronic shell structure may be 
determined from the spectrum of a one-particle Hamil- 
tonian, provided the V(r) is. chosen sufficiently close to 
the self-consistent potential.0 

For small sodium clusters of 8 up to 92 atoms it has 
been shown that the shell structure arises from succes- 
sively filling the 2(21 + l)-fold degenerate-.energy levels 
£n : i of a spherically symmetric potential.E3 To uniquely 
characterize the energy levels for such potentials we use 
the number n of nodes in the radial wave function, while I 
denotes the angular momentum quantum number x3 Go- 
ing to larger clusters reduces the influence of a single 
energy level on the shell structure. Nevertheless even for 
clusters of several thousand-atoms electronic shell struc- 
ture has been observed. This suggests that the en- 
ergy levels bunch together in groups of (nearly) constant 
energy. Actually, the experimental results-, for sodium 
clusters of up to approximately 1500 atomstTllJ indicate 
that these bunches may be formed by degenerate energy 
levels e ni i corresponding to definite values of 3n+Z. It has 
been pointed out that quantities of the type an + 1 char- 
acterize two exactly solvable quantum mechanical prob- 
lems: the hydrogen atom, in which levels with the same 
n + / are degenerate, and the isotropic harmonic oscilla- 
tor, where the degeneracy is described by the quantum 
number 2nJ~L— With this analogy in mind it has been 



conjecture that there might be a set of solutions of 

the radial Schrodinger equation displaying the degener- 
acy described by the quantum number 3n + I. 

At first sight a theorem from classical mechanics seems 



to imply that there is no such radial potential. The the- 
orem states that only for central potentials_of the form 
— 1/r and r 2 all bounded orbits are closed£3 Assuming 
that this property is caused by a dynamical symmetry of 
the system which uniquely corresponds to a symmetry of 
the quantized system, this would imply that only these 
two problems exhibit "accidental" degeneracies. Unfor- 
tunately the situation is not that simple. While it is 
known that the degeneracies for the quantum mechani- 
cal Kepler problem and the harmonic oscillator are a con- 
sequence of the dynamical symmetries so (4) and su(3), 
respectively, it turns out that all classical systems involv- 
ing 3 degrees of freedom automatically possess both so (4) 
and su(3) symmetry.E3 To our knowledge, the question 
whether the hydrogen atom and the harmonic oscilla- 
tor are the only radially symmetric quantum mechanical 
systems possessing a dynamical symmetry has not been 
settled so far. A discussion of near degeneracies of the 
type an + (31 in the context of semiclassiG|al quantization 
has been given by Bohr and Mottelson.E2l 

Besides being interesting in its own right, identifying 
an exactly solvable problem possessing a quantum num- 
ber 3n + I would be important for two reasons. First, 
it would provide us with the optimal basis set for self- 
consistent computations of the total energy of metal clus- 
ters. Second, it would give the ideal starting point for 
treating realistic cluster potentials by perturbation the- 
ory. 

In the present paper, we investigate whether the num- 
ber 3n + I is analogous to the quantum numbers n + I 
and 2n + I. First we argue that a degeneracy described 
by the quantum number 3n + I cannot be caused by a 
dynamical symmetry related to a classical Lie algebra. 
Next we show that in three dimensions there are exactly 
two potentials fulfilling the shape-invariance condition of 
supersymmetric quantum mechanics for spherically sym- 
metric problems, namely the Coulomb potential and the 
harmonic oscillator. Finally we numerically optimize ra- 
dial potentials to exhibit the 3n + /-degeneracy. We find, 
within the limits of our calculations, that there is no ex- 
actly solvable 'cluster potential' with 3n + /-degeneracy. 
We therefore conclude that 3n+l is not a proper quantum 
number as n + I and 2n + I are. 
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system 


quantum number 


symmetry 


hydrogen atom 


n + I 


so (4) 


harmonic oscillator 


2n + l 


su(3) 


cluster 


3n + l 


? 



TABLE I. Possible analogy between systems with a dy- 
namical symmetry and a conjectured 'cluster problem', being 
characterized by a quantum number 3n + I. 



II. CLASSICAL SYMMETRY 

A quantum mechanical system is said to exhibit a sym- 
metry if its Hamiltonian is invariant with respect to the 
operations of some Lie algebra. A well-known example 
is spherical symmetry: the Hamiltonian H of a radially 
symmetric system commutes with the angular momen- 
tum operators: 



[L i ,H} = 0. 



(1) 



Since {L x , L y ,L z } forms a Lie algebra, namely the angu- 
lar momentum algebra su(2), the Hamiltonian is said to 
be invariant under su(2). As an immediate consequence 
of (0) the eigenspaces 7i €ji , of the Hamiltonian are in- 
variant under the action of the Li, i.e. they are repre- 
sentation spaces of the symmetry algebra. For a generic 
radial potential, the Ti. en , will be irreducible representa- 
tion spaces of the angular momentum algebra. For these 
systems su(2) is the full symmetry. 

There are, however, special cases in which the eigen- 
spaces of the Hamiltonian are reducible. Examples are 
the hydrogen atom and the isotropic harmonic oscillator. 
For the former, energy levels e n j having the quantum 
number n + / in common are degenerate, for the latter 
the respective quantum number is 2n + l. It is well known 
that these degeneracies are not accidental, rather they 
are the consequence of a symmetry higher than spheri- 
cal. For the hydrogen atom this symmetry is described 
by the algebra so (4) while for the harmonic oscillator the 
symmetry algebra is the su(3).lj Having identified 3n + l 
as a candidate for a new quantum number, the question 
arises as to whether the corresponding degeneracy is also 
due to some hidden symmetry. 

To answer this question, we exploit the fact that 
eigenspaces of a Hamiltonian are representation spaces 
of its symmetry algebra. In particular, the dimensions of 
these s pac es must coincide. Fortunately, representation 
theory! 22 !! 23 ! provides us with all the information needed. 
The irreducible representations of the classical Lie al- 
gebras (su(/c + 1), so(2fc + 1), sp(2k), and so(2/c), cor- 
responding to the algebras Ak, B}~, Ck, and Dk in the 
complete classification scheme of all semisimple Lie alge- 
bras) can be labeled using the Cartan subalgebra. For a 
Lie algebra of rank k this is a set of k commuting oper- 
ators. They can be chosen such that they have integer 



eigenvalues. For any given representation one can find a 
basis which diagonalizes the operators in the Cartan sub- 
algebra. Each basis vector is then characterized by the 
corresponding integer vector of eigenvalues (mi, . . . , m^), 
its 'weight'. Furthermore each irreducible representation 
is uniquely determined by the highest weight of its ba- 
sis vectors according to lexicographical ordering. The 
dimension of the irreducible representation with highest 
weight (m, . . . , nk) is given by Weyl's dimensionality for- 
mula 



dim 



(r,c(ni, . 



n 

aeA 4 



(<*) 

3 = 1 



k 



,(«) 



(2) 



The constants c^' and the set A + of positive roots only 
depend on the type of the Lie algebra, not on any par- 
ticular representation. Hence for a given classical Lie al- 
gebra, the dimensions of the irreducible representations 
are given by a polynomial in the highest weights. As 
an example we write down explicit expressions for the 
dimensions of the classical algebras of rank two: 



dim FT A 2 (ni,n 2 )j = ■ 

dim(r B2 (ni, n 2 )j 

(ni + l)(m 



(ni + l)(ni+n 2 + 2)(n 2 + l) 



(3) 
(4) 



n 2 + 2)(2ni+n2 + 3)(n 2 + l) 



G 



dim(r D2 (ni,n 2 )j = (ni + l)(ri2 + 1). 



(5) 



These have to match the dimensions of the eigenspaces of 
the Hamiltonian. Given the quantum number N associ- 
ated with the degeneracy of energy levels e n ,i, the orders 
gN of the degeneracies can be readily computed. For the 
quantum number N\ = n + I we obtain 



N 2 = 2n 



9N ± 

I leads to 



(JVi 



1) 



9n 2 - 

and N$ = 3n + I to 



(N 2 + l)(N 2 + 2) 



(G) 



(7) 



9N 3 



i(7V 3 + 1)(JV 3 



(8) 



for N3 mod 3 = 1, 
3), otherwise. 

Through comparing with expressions (§)-(|^), we find 

g Nl =dim(r i52 (iV 1 ,iV 1 )), (9) 

g N2 =dim(r^(AT 2 ,0)), (10) 



indicating that the eigenspaces of the Hamiltonian for 
the hydrogen atom might be irreducible representation 
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spaces of the Lie algebra D% = so (4) with highest weights 
(711,712) = (N±,Ni) — as indeed it is. Similarly for 
the harmonic oscillator the eigenspaces could be the irre- 
ducible representation spaces of A2 = su(3) with highest 
weights (711,712) = (iV2,0) — again, as indeed it is. For 
the quantum number 3n + I the situation is fundamen- 
tally different. The degeneracies of the eigenspaces are no 
longer given by a closed formula, i.e. there is no simple 
relation between the quantum number and the highest 
weights of some Lie algebra. If there existed a quantum 
mechanical operator corresponding to the quantum num- 
ber, which would be a Casimir operator of the symmetry 
algebra, then such a simple relation should exist, since 
the eigenvalues of a Casimir operatot-are given by ratio- 
nal functions in the highest weights.EJ Thus, we cannot 
find a classical Lie algebra that could cause a degeneracy 
corresponding to the quantum number 3ti + 1. The same 
reasoning holds for quantum numbers an + 1 with a > 3. 



III. SUPERSYMMETRY 

An elegant approach to exactly solvable problems, is 
provided by supersymmetric quantum mechanics Hj'Ea 
The basic idea of the method is that a given Hamilto- 
nian 



H' 



d 2 
dx 2 



(11) 



can be factorized if the energy of the ground state (x) 
is set to zero: 



H~ = Q+Q- , where Q 



± _ 



^ dx 



(*o)' 

*r7 



(12) 



Here the prime denotes the derivative with respect to x. 
Reversing the order in the factorization, a new Hamilto- 
nian H + = Q~Q + can be defined. It is called the super- 
symmetric partner of H~. Obviously the intertwine 
the Hamiltonians: 



H-Q+ = Q + H A 
H+Q- = Q-H- 



(13) 
(14) 



Thus, except for the ground state of H~ , which is anni- 
hilated by Q~, the spectra of the Hamiltonians H T are 
identical. This situation is illustrated in Fig. [I]. 

The potentials corresponding to the Hamiltonians H T 
can be derived from the superpotential 



W(x) = -<M. 



V T (x) = (\V(x) 2 =F W'(x) 



(15) 



(16) 



For a generic Hamiltonian the potential V (x) and its 
supersymmetric partner V + {x) will be of different shape. 



H- 



Q+Q- 



H+ = Q-Q+ 




FIG. 1. Spectra of supersymmetric partner Hamiltonians 
H~ and H + (non broken supersymmetry) . 

If, however, they only differ in some parameters and a 
constant offset, i.e. if 



V + (a ;x) = V-{a 1 ,x)+R(a ), 



(17) 



the corresponding Hamiltonian is exactly solvable.0 In 
what follows we shall refer to the expression in equation 
( |l7| ) as the shape-invariance condition. In such a case, 
starting from the potential V~(ao',x), a whole sequence 

V-(a k ;x) ^V+ia^x) = V~ (a k+1 ;x) + R(a k ) (18) 

of supersymmetric partner potentials can be constructed. 
From the knowledge of the ground-state energy in the 
potential V~(a; x) as a function of the parameters a, the 
whole spectrum of the initial Hamiltonian H~ can be 
obtained: 



n-l 



E- = Y J R{ak). 



(19) 



fc=0 



We now investigate whether a radial problem with de- 
generacies according to the quantum number 3n + / can 
be constructed using the above ideas. More specifically, 
we are looking for any potential U(r) such that energy 
levels with different angular momentum / are degenerate. 
For the radial potentials 



V{l;r) = U(r) + 



1(1 + 1) 



the shape-invariance condition ( |l7| ) reads 

V + (l;r)=V-{l;r)+R(l). 
Using ( p0| ) this can be rewritten as 

V+(l;r) = V-(l;r) + ^- + 20i, 

where we have introduced the notations 

2Xi = 1(1+1) -1(1 + 1) and 2/3i = R(l). 



(20) 



(21) 



(22) 



(23) 



3 



Inserting ( |l6| ) into (p2|), we obtain a differential equa- 
tion for the superpotential W(l\r), the general solution 
of which is 

W(l;r)=-^+ai+pir. (24) 

The corresponding potentials are 

Ai(AiTl) 2X iai 



(25) 



-(2A / ±l)/3 / + (a / +/3 / r) 5 



The choices A/ = 1 + 1 and A; = — I will make V _ (Z;r) 
a radial potential with angular momentum I. To ensure 
that we can normalize the ground state, the parameters 
in the above equations must satisfy 



A; > — — , and 



Pi > 0, 
or 

Pi = 0, ai > 0. 



(26) 



We thus have to choose \i = I + 1 > — 1/2. Inserting this 
into equation (|25|) yields 



V + (l;r)-V-(l + l,r) (27) 
= {(a; + Pirf - (aj+i + A+ir) 2 } 

'2(J + l)ai 2(Z + 2)a /+1 



+ { - (2/ + l)A + (2/ + 5)A+i}- 



To be consistent with ( |l7| ) , we require the right-hand side 
of the above equation to be constant. By comparing the 
coefficients, we find exactly two solutions of the shape- 
invariance condition: 



ai+i 
ai 



l + l 
1 + 2 



and pi = 0, 



on = and Pi 



const. 



(28) 



(29) 



The first solution corresponds to the Hamiltonian for the 
hydrogen atom. Setting ai = q/{2(l + 1)) we obtain the 
energy spectrum 



e n;l = E(l)- = 



(30) 



4(n + 1 + l) 2 

pertaining to the Coulomb potential V(r) = —q/r. Writ- 
ing the constant in the second solution as Pi = uj/2 
we find the spectrum of the harmonic oscillator V(r) = 

e„,z = E(l)~ = (2n + Z + 3/2)u/. (31) 

Obviously, the spectra of these systems are characterized 
by degeneracies according to the quantum numbers n + I 
and 2n + l. Thus we have shown that there exist exactly 
two physically acceptable radial potentials which fulfill 
the shape-invariance condition (21), namely those asso- 
ciated with the quantum numbers n+l and 2n+l. It thus 
appears that these two quantum numbers are unique. 



IV. NUMERICAL OPTIMIZATION 



Clearly, arguments of the type given above can only 
exclude the existence of a system with quantum number 
3n + I under the assumption of some specific mechanism 
(e.g., classical symmetry or supersymmetry) that would 
cause the degeneracy. Thus, by such reasoning it is not 
possible to give a general proof that such a system can- 
not exist. Furthermore, it is conceivable that 3n + I is 
only an approximative quantum number, i.e. that levels 
of quantum number 3n + I are not exactly but merely 
near-degenerate. Therefore, we will now approach the 
problem at hand from a completely different direction. 

We try to construct cluster-potentials exhibiting de- 
generacies according to 3n + 1 by numerical optimization. 
The basic idea is to transform the continuous eigenvalue 
problem for the Hamiltonian 



H 



2m 



A + V(r) 



(32) 



into a finite dimensional problem by discretizing the 
Schrodinger equation. Discretizing the potential 



{V(r) | r E [0,oo)} h-> {v 4 | V, = V(n)} 



(33) 



the eigengenvalues can be considered as functions of a 
finite number of potential parameters Vi instead of being 
junctionals of the potential: 



V{r) 



(34) 



A potential with the postulated in + I degeneracy can 
then be found by numerically optimizing a suitable initial 
potential with respect to a cost function, which measures 
the deviation from the desired degeneracy. By a straight- 
forward application of the above strategy we have found 
that the V^s for successive i's oscillate wildly. For such 
potentials {Vi} the approximation (|34|) is clearly invalid. 
Thus we have imposed a constraint on the V^'s to keep 
the potential smooth. Furthermore we have made sure 
that the number of electrons the potential can hold for a 
given Fermi energy does not vary too much. 

We have performed such optimizations starting from 
Woods-Saxon potentials 



V(r) 



-Vo 



exp(^) 



(35) 



with sodium-like parameters (Vo = 0.46 Ry and a = 
0.94 ao) for varying cluster sizes Rq. For the discretiza- 
tion we have used a mesh of step width Ar w r s /40. 
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FIG. 2. Result of an optimization run for the quantum 
number 2n + I. The optimized potential (full line) along 
with the energy levels (horizontal lines) are plotted. The 
computation was started from a Woods-Saxon potential of 
width _Ro = 8r s , with r s = 3.93 ao (dotted line). The 
quadratic potential (dashed-dotted line) exhibiting the exact 
2n + /-degeneracy is given for comparison. 

A cost function has measured the non-degeneracy of the 
energy levels, the smoothness of the potential and the 
variation in the number of electrons that fit into the 
potential for a given Fermi energy. To find the non- 
degeneracy we have determined for each quantum num- 
ber N3 = 3n + I the average energy 



e N 3 



E 



3n+l=N 3 



2(2/ + 1) e n ,i 



E 



3rM 



= JV; 



2(2/ + 1) 



(36) 



The deviation of the eigenstates from the desired degen- 
eracy was then defined as 



E 

N 3 



£JV 3 



E 



3n+l=N 3 



2(2/+ 1) 



(37) 



As a measure of the roughness of the potential repre- 
sented by the Vi we have used the expression 



(38) 



which disfavors sharp bends in V(r). The number of 
electrons fitting into the potential was estimated by the 
Thomas-Fermi expression 



3-7T 



3/2 



(39) 



where r out is the classical turning point for the potential 
V(r). 

Since the Woods-Saxon potential has only a finite num- 
ber of bound states, the summation ([37]) has to end at 
some maximal N, which we will denote by N max . 

To demonstrate that with the above measures we in- 
deed can find potentials with a given degeneracy, we have 
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FIG. 3. Potentials for increasing cluster radius 7?o, op- 
timized for degeneracies according to the quantum number 
3n + /. The optimization was started from Woods-Saxon 
potentials with sodium-like parameters. r 3 = 3.93 ao is the 
Wigner-Seitz radius of Na. The lines indicate the position of 
the energy levels. 
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FIG. 4. Optimization of cluster potentials for degeneracies according to the quantum number 3n + l. On the left, potentials 
are plotted along with the energy levels. The grouping of levels is shown on the right, where the eigenenergies t n ,i are plotted 
over the corresponding angular momenta I. Levels having the same quantum number 3n + l are connected by a dotted line. The 
plots at the top show the spectrum of the Woods-Saxon potential, which has been used as initial potential for the optimization 
runs shown in the plots below. These runs only differ in the number N max of shells included in the cost function. 



tested our program on the quantum number 2n + I. We 
start from the sodium-like Woods-Saxon potential that 
we also use for the 3n + /-optimizations. Although this 
initial potential is much closer to a 3n + I than to a 
2n + /-degeneracy, it turns out that during optimization 
the energy levels are very easily arranged into equidistant 
shells. For all choices of N max nearly perfect degeneracy 
is achieved for the shells with N < N max . Interestingly, 
even the levels which are excluded from the cost function 
(i.e. those above N max ) tend to form degenerate shells. 
Figure || shows the energy spectrum along with the po- 
tential resulting from the optimization of a Woods-Saxon 
potential with Rq — 8r s . As can be seen the optimized 
potential approximates well the quadratic potential (in- 
dicated by the dash-dotted line) associated with the 2n+l 
degeneracy, although it-is very different from the initial 
potential (dotted line).E3 

The results of the optimizations for the quantum num- 
ber 3n + l are quite different. They are shown in figure 0. 



To make the results for different cluster sizes comparable, 
we have chosen N rnax to correspond always to the shell 
closest to Fermi energy. It turns out that for small clus- 
ters the energy levels can be forced to form nearly degen- 
erate shells according to the quantum number 3n+l. The 
corresponding potentials, however, exhibit large oscilla- 
tions inside the cluster. These seem to be fairly unrealis- 
tic, when compared to those obtained form self-consistent 
calculations For larger clusters the situation becomes 
even worse. Here, although the optimized potentials os- 
cillate strongly, degeneracy can merely be achieved for 
the states with in + / just below N max . In other words, 
only just below Ef a grouping of almost degenerate en- 
ergy levels can be observed. Above and further below 
N m ax there are no large gaps in the spectrum, and hence 
no shells. 

The above results suggest that the outcome of a de- 
generacy optimization strongly depends on N rnax , indi- 
cating that there is no unique 'cluster potential' which 
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could be approximated by the optimized potentials. To 
demonstrate this dependence, we have performed calcu- 
lations with different values of N max for a fixed cluster 
size. The results are shown in Fig. |J. At the top the 
Woods-Saxon potential, which was used as starting point 
in the optimization runs, along with its energy spectrum 
is plotted. On the left-hand side the potential is shown 
as in Fig. [| (except that we restrict the energy axis to 
a smaller range). To the right the eigenstates e n j are 
arranged according to their angular momenta States 
belonging to the same shell -/V3 = 3n + I are connected 
by a dotted curve. A shell N3 is (near-)degenerate if this 
curve is a horizontal line. The greater its slope, the worse 
the degeneracy. Optimizing the potential for two differ- 
ent values of N max results in the potentials shown in the 
lower part of Fig. [|. These potentials are markedly differ- 
ent. Not only the amplitude of the oscillations, even the 
number of maxima in V(r) changes with N max . As can 
be seen from the spectra on the right hand side the dotted 
curves connecting states, which belong to the same shell 
N3 = 3n + I are concave for N3 > N max and convex for 
N3 < N max . Comparing the spectra for N max = 10 and 
12, we observe that there is a tradeoff between straight- 
ening out the dotted lines above N3 = 10 and bending 
those corresponding to smaller quantum numbers; i.e. it 
seems to be impossible to obtain the desired degeneracy 
for all shells. Furthermore it turns out that the higher 
N m ax (he. the more states are included in the cost func- 
tion), the stronger are the oscillations in V(r). 

If a 'cluster potential' having a shape similar to a 
Woods-Saxon potential existed, then the optimized po- 
tentials would have approximated that unique poten- 
tial. Like in the optimizations for the quantum number 
2n + I, the degeneracies would become the better the 
larger N max . This is clearly not what we have found 
in our calculations. Instead, the optimized potentials 
strongly depend on N max and further we have not found 
any potential whose complete spectrum is degenerate ac- 
cording to the quantum number 3n + I. 



V. CONCLUSIONS 

By use of analytic methods as well as numerical calcu- 
lations we have argued that there is no 'cluster problem' 
analogous to the exactly solvable problems of the hydro- 
gen atom and the harmonic oscillator. The sequence of 
quantum numbers n + / and 2n + 1 characterizing the lat- 
ter systems thus cannot be simply extended to quantum 
numbers 3n + 1. In this sense it seems that the quantum 
numbers n + I and 2n + I are unique. 

In spite of this, the higher 'quantum numbers' an + (31 
do have some meaning. The analysis of the energy lev- 
els near the Fermi energy found in self-consistent jellium 
calculations reveals that the electronic shell and super- 
shell structure mainly arises from an intesplay of near 
degeneracies of the type 3n + 1 and An + 1 SB As has been 



pointed out by Bohr and Mottelson,E2l such near degen- 
eracies of energy levels corresponding to numbers an + (3l 
are related to classical periodic orbits. For example, tri- 
angular orbits are connected with 3n+Z-degeneracies, and 
square orbits with An + 1. This simple correspondence is 
however lost in a full semiclassical analysis. Nevertheless 
for typical cluster potentials a periodic orbit expansion!!! 
shows that the electronic shell and supershell structure is 
essentially captured by only considering the most impor- 
tant contributions which stem from the triangular and 
the square orbits. Using the simple picture of Bohr and 
Mottelson this is consistent with the conjectured degen- 
eracies as well as with those that can be found in jellium 
calculations. 

Thus, the structures observed in the mass spectra of 
warmed metal clusters are intimately related with the 
quantities 3n + I and An + I. These numbers are, how- 
ever, not analogous to the quantum numbers known from 
the exactly solvable models. 
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